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Abstract 

In this paper, we study dimension reduction of the three-dimensional (3D) 
Gross-Pitaevskii equation (GPE) modelling Bose-Einstein condensation un- 
der different limiting interaction and trapping frequencies parameter regimes. 
Convergence rates for the dimension reduction of 3D ground state and dy- 
namics of the GPE in the case of disk-shaped condensation and cigar-shaped 
condensation are reported based on our asymptotic and numerical results. In 
addition, the parameter regimes in which the 3D GPE cannot be reduced to 
lower dimensions are identified. 

1 Introduction 

Since its realization in dilute bosonic atomic gases [1, 16, 11], Bose-Einstein con- 
densation (BEC) has been produced and studied extensively in the laboratory 
[27, 29, 20], and has afforded an intriguing glimpse into the macroscopic quantum 
world [29]. The experimental advances in BEC have spurred great excitement in 
the atomic physics community and renewed the interest in studying the collective 
dynamics of macroscopic ensembles of atoms occupying the same one-particle quan- 
tum state. Needless to say that this dramatic progress on the experimental front 
has stimulated a wave of activity on both the theoretical and the numerical front. 

At temperatures T much smaller than the critical temperature T c [22] , a weakly 
interacting BEC is well described by the macroscopic wave function if> = ip(x,t) 
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whose evolution is governed by a self-consistent, mean field nonlinear Schrodinger 
equation (NLSE) known as the Gross-Pitaevskii equation (GPE) [19, 28, 29] 



ih 



dip(x,t) h 



dt 



_ — V>(x, t) + V(x)^(x, t) + NUqMx, t) |>(x, t), (1.1) 
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where x = (x, y, z) is the spatial coordinate, m is the atomic mass, h is the 
Planck constant, N is the number of particles in the condensate, V(x) is an ex- 
ternal trapping potential. When a harmonic trap potential is considered, V(x) = 
Y (oo 2 x 2 + oo 2 y 2 + uJ 2 z 2 ^j with u x , oo y and oo z being the trap frequencies in x, y and 
^-direction, respectively. Without loss of generality, we assume uj x < uj y < uj z . 
jj _ Anh a s dggcribgg the interaction between atoms in the condensate with a, the 
s-wave scattering length. The wave function is normalized according to 

|^(x,t)| 2 dx=l. (1.2) 

Following the physical literatures [12, 7, 5, 29], introducing the rescaling: t — > 
x — > aox with ao = \/^~, an d ip — > we get the following dimensionless GPE 



under the normalization (1.2) in 3D: 



where /3 



C/qAT 
aphDx 



~VV(x, t) + V(x)^(x, t)+(3 |V(x, t) | 2 ^(x, t), (1.3) 
and V(x) = § (x 2 + 7 2 y 2 + 7 2 ^ 2 ) with % = ^ and 
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Iz = 



To find the stationary solution of (1.3), we write 

^(x,t) = 0(x) e"**, 



(1.4) 



where /x is the chemical potential of the condensate and 0(x) is a function indepen- 
dent of time. Substituting (1.4) into (1.3) gives the following equation for (/j, 0(x)): 



/i <Hx) = --V 2 0(x) + V(x)0(x) + /5|0(x)|V(x), 



x G 



(1.5) 



under the normalization condition (1.2) with ip = <fi. 

This is a nonlinear eigenvalue problem with a constraint and any eigenvalue \i 
can be computed from its corresponding eigenfunction 0(x) by 



where 



= /i(0) = I 
Jr- 



: |V0(x)| 2 + \/(x)|0(x)| 2 + /3|0(x) 



dx = £(0) + £ int (0), (1.6) 



|V0(x)| 2 + \/(x)|0(x)| 2 + f|0(x)| 4 



(1.7) 



E ^ = L 

Jr- 

Emt(4) = ^f I0(x)| 4 rfx, E kin (<f>)=[ i|V0(x)| 2 rfx. (1.8) 
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The ground state of a BEC is usually defined as the minimizer of the following 
minimization problem: 

Find (fj, g , 4> g E S) such that 

E g := E{4> 9 ) = min£(0), fi g := (ji(<f> g ) = E(<f> 9 ) + E int (<f> g ), (1.9) 

where S = {0 | ||0|| = 1, E{<p) < oo} is the unit sphere. The existence of a unique 
positive minimizer of the minimization problem (1.9) is given in [25]. 

In an experimental setup, the trapping frequencies in different directions can be 
very different. Especially, disk-shaped and cigar-shaped condensation were observed 
in experiments [9, 15, 24, 29]. The 3D GPE (1.3) is formally reduced to 2D GPE 
in disk-shaped condensation and to ID GPE in cigar-shaped condensation in the 
literatures [15, 26, 31, 24, 29, 25, 23, 7, 5, 3, 30]. Mathematical and numerical 
justification for the dimension reduction of 3D GPE is only available in the weakly 
interaction regime, i.e. (3 = o(l) [6, 10]. Unfortunately, in the intermediate or 
strong interaction regime, no results are available. The aim of this paper is to study 
dimension reduction of the 3D GPE (1.3) under different limiting parameter regimes. 
In addition, we will provide convergence rate for the dimension reduction based on 
our asymptotic and numerical results. 

The paper is organized as follows. In section 2 we study dimension reduction for 
the ground state of the 3D GPE (1.5) with replusive interaction (3 > and provide 
convergence rates. Convergence rate for dimension reduction of 3D time-dependent 
GPE (1.3) is reported in section 3. Finally some conclusions are drawn in section 4. 



2 Dimension reduction for ground states 

In this section, we will discuss dimension reduction for BEC ground states and derive 
approximate ground states as well as their energy and chemical potential for the 3D 
GPE (1.5) under different parameter regimes of (3, ^ y and j z . 

2.1 Isotropic shaped condensation 

In the case of isotropic condensates, i.e. = 0(1) and j z = 0(1) (■<=>• co y ~ u x 
and uj z ~ u x ), there are three typical regimes: 

Regime I. Weakly interacting regime, i.e. (3 = o(l), the ground state is well approx- 
imated by the harmonic oscillator ground state [13, 15, 29, 7]: 

g (x) ^0 ho (x,y,z) = ^y^l'\ - x2+lv i +1 ^\ xGR 3 , (2.1) 
E g » £(0 ho ) = ± (1 + ly + lz ) + 0((3), \0\ « 1, (2.2) 

Hg » M0ho) = \ (1 + ly + lz) + O(P). (2.3) 
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Regime II. Intermediate interacting regime, (5 = 0(1), the ground state can be 
obtained by solving the 3D minimization problem (1.9). Different numerical methods 
were proposed in the literatures for computing the ground states [7, 4, 3, 12, 13]. 

Regime III. Strong interacting regime, (5 >> 1, noticing (A. 7) and (A. 6), the ground 
state is approximated by the Thomas- Fermi (TF) approximation [13, 15, 29, 7]: 



09( x) „ c(x) = ( \[W F - /ft n*) < 

[ 0, otherwise, 



(2.4) 

otherwise, 

C = i(^) 2/5 - < 2 ' 5 > 

» 7< + |js 0» + G.) = ^ F + O (^) , (2.6) 

ft » c + ^? cm + G »> = "» TF + ° ■ < 2 ' 7 > 

For 7 y = 7 2 = 1, (2.6) and (2.7) were confirmed numerically in [7]. 
2.2 Disk-shaped condensation 
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In the case of disk shaped condensates, i.e. ^ y = 0(1) and ^ z 3> 1 (< 
and o; 2 ^> u; x ), we set 

7z ^y 1 / 4 7 2 z 2 

/^~/^ + y, 9 (x) ps 0(x,y)0 ho (z) with ho (z) = 2 . (2.8) 

Plugging (2.8) into (1.5), multiplying both sides by (f)ho(z) and integrating over 
z G (— oo, oo), we get 

H<j>{x,y) = ~A(j>(x,y) + V 2 {x,y)<j){x,y) + ft \4>{x, y)\ 2 <f)(x, y), (x,y) el 2 , 

(2.9) 

where 

V 2 (x, y) = \ (x 2 + rfy 2 ) , ft = (3 r |0 ho (^)| 4 dz = (3J^. 

Using the results in Appendix A for 2D GPE, again we get approximate ground 
state in three typical regimes: 



Regime I. Weakly interacting regime, i.e. ft = (3y r y z /27t = o(l), the ground state 
is approximated by the harmonic oscillator ground state: 

9 (x) ps (f) ho (x, y)(j) ho (z) = (j) ho (x, y,z), x G R 3 , (2.10) 

E g ps | + l±2k + o(ft) = | + i±2k + o (/3 72 1/2 ) , (2.11) 
^ - y + ^ + 0{ft) = I + + O (HJ 2 ) , 7, » 1 & & a = o(l). (2.12) 
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i 1 z 




1/25 


1/100 


1/400 


1/1600 


8 = 


1 


2 23E-3 


9.97E-4 


4.17E-4 


1.68E-4 


rate 






0.58 


0.63 


0.65 


8 — 


10 


1.23E-2 


4.72E-3 


1 74E-3 


5 80E-4 


rate 






0.70 


0.72 


0.80 


P = 


100 


4.46E-2 


1.64E-2 


5.93E-3 


2.13E-3 


rate 






0.72 


0.73 


0.74 


P = 


1000 


1.38E-1 


5.31E-2 


1.94E-2 


6.99E-3 


rate 






0.69 


0.73 


0.74 


P = 


10000 


3.46E-1 


1.58E-1 


6.19E-2 


2.27E-2 


rate 






0.56 


0.68 


0.72 



Table 1: Error analysis of \\4> g (-) — 0g S (-)IU 2 f° r the ground state in 3D with a 
disk-shaped trap. 



Regime II. Intermediate interacting regime, i.e. = 0(1), the ground state can be 
approximated by 

(2.13) 



9 (x) » 0f (x) := 4>f{x,y)^ {z), x G M 3 . 



;2D 



E n E 



DS 



/i 9 W /i 
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DS 



= E^(x,y)(f> bo (z)) 
/ x(0f(a:, 2 /)0 ho (z)) = 



| + ^ 2D (0f) 
T + ^D(0f) := 



- + #„ 2D , 



— + « 2D 



(2.14) 
(2.15) 



where 



E; 



2D 



|V0f| 2 + y 2 (x,y)|0f| 2 



+ fl0fl 4 



(.2D 1 4 



dxdy, 



dxdy. 



Here 2D , _E| D and /i 2D are the ground state, energy and chemical potential of the 
2D problem (2.9). In this case, one only needs to solve a 2D problem numerically 
and thus significantly save computational time, memory and cost. 

To verify (2.13), (2.14) and (2.15) numerically, Table 1 list the error ||0 9 (x) — 
DS (x)|| £ 2, and Figure 1 shows the error and \E g — E^ s \, with j y = 1 for different 
(5 and j z . Here and in the following, the ground state <p g is computed numerically 
by the continuous normalized gradient flow (CNGF) with a backward Euler finite 
difference (BEFD) discretization [4]. 

From Tab. 1, Fig. 1 and additional numerical results [18], when ^ z 1, (3% — 
P^lz/^n = 0(1) or 1, and j3% 3 ^ 2 = o(l), we can draw the following conclusion: 



110,- 



iDSl 



i 2 



O 



3/4 

1* 



DS\2i 



110 - (O 



L 1 



o 



'(3 1 ' 2 In 
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On the contrary, when 7 Z > 1, = I3^ z /2ti = 0(1) or > 1, and /?77 3/2 = 0(1) 
or ^> 1, the errors do not decrease when 7 Z increases. This suggests that one 
cannot reduce the 3D GPE (1.5) to 2D GPE (2.9) when (3 > 1 and ^ z > 1 with 



Regime III. Strong interacting regime, (3% = flylz/Zn ^ 1, noticing (A. 17), the 
ground state is approximated by the multiplication of the TF approximation in 
xy-plane and the harmonic oscillator approximation in z-direction: 



9 (x) » 9 1M (x) := ^(x,y)0 ho (z), x G (2.16) 

where 

aTF (x )= i yf(fJ$E - V 2 (x,y))/P$, V 2 (x,y) < //^, (2 U) 

[ otherwise, 

foa \V2 ol/2 1/2 1/4 
^ = \^— J - 2 V%3/4 • (2- 18 ) 

Plugging (2.13), (2.9), (A.7) with d = 2, (2.18), (A.16) with d = 2 and & = (3% = 
P\f^ (2-14), we get the approximate energy 

E g = E( ( f )g ) = E(<p?(x,y)<PUz)) + o(^^ 
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1/7* 


1/25 


1/100 


1/400 


1/1600 


8 = 1 


5.45E-1 


4.44E-1 


3.57E-1 


2.88E-1 


rate 




0.15 


0.16 


0.15 


p = 10 


2.66E-1 


2.15E-1 


1.74E-1 


1.40E-1 


rate 




0.15 


0.16 


0.15 


/3 = 100 


1.29E-1 


1.03E-1 


8.43E-2 


6.77E-2 


rate 




0.16 


0.14 


0.16 


/3 = 1000 


1.40E-1 


6.49E-2 


4.12E-2 


3.19E-2 


rate 




0.55 


0.33 


0.18 



Table 2: Error analysis of \\(f> g — TF1 1| z, 2 for the ground state in 3D with a disk-shaped 
trap. 



= £ + s D (0 + °'" ln7 - 



1/2 
1* 



2 3tt 3 / 4 + 2(/3 2 72 ) 1 /4 



ln(/3 2 7z ) +2G 2 -ln27r 



+ 



7 



1/2 



= £j F1 + 



/ln(/3 2 7z ) , /31n 7 



+ 



V(/? 2 7,) 1/4 7 , 1/2 



C2-19) 



where 



TTiTFl I* 



lz 



2 3/ V 2 (l 2 7,) 1/4 
3tt 3 / 4 

Similarly, we get the approximate chemical potential: 

ln((3 2 lz ) , /31n 72 



A*9 ~ A* s 



TF1 + 



.(/5 2 7,) 1/4 7 i /2 



where 



TF1 = Jl + 
2 + 



il l2 ^ 2 izY' A 



^2.20) 



(2-21) 



^2.22) 



"» 2 ' 2^4^3/4 • 

To verify (2.16), (2.19) and (2.21) numerically, Tables 2 and 3 list the errors 
||0 9 (x) — TF1 (x)|| i 2 and ||0 9 (x) — TF1 (x)||ioo, respectively, and Figure 2 shows the 
error \E g — -Ej F1 |, with 7j , = 1 for different /3 and 7z . 

From Tabs. 2&3, Fig. 2 and additional numerical results in [18], when 7y = 
0(1), 7z > 1, /3£ = (3^ z /2tt > 1 and /3 7 ~ 3/2 = o(l), we can draw the following 



conclusion: 



114 



Ml - (0 TF1 )^ 



IL 1 



Q ( CiP) ln 7z 



1/4 

7* 



1/7* 1/25 1/100 1/400 1/1600 



(3 = 1 0.3880 0.3148 0.3107 0.3749 
13 = 10 3.845E-2 4.589E-2 5.337E-2 6.198E-2 
(3 = 100 6.368E-3 7.048E-3 8.597E-3 9.757E-3 



Table 3: Error analysis of \\(f> g — tf1 ||l°° for the ground state in 3D with a disk- 
shaped trap. 




Figure 2: Convergence rate of \E g — -Ej F1 | in 3D with a disk-shaped trap: (a) With 
respect to 7*; (b) with respect to (3. 



where C{(3) depends on (3. These results confirm the asymptotic results (2.19) and 
(2.21). Furthermore, the additional numerical results in [7, 18] indicate that TF1 (x) 
does not converges pointwise to the ground state 9 (x) when j z — > 00 and (3 > 
(cf. Tab. 3). 



2.3 Cigar-shaped condensation 

In the case of cigar shaped condensates, i.e. 7 y ^> 1 and j z ^> 1 (■<=>• u> y ^> uj x and 

uj z ^> u> x ), we set 



vl/4 



V o * . 0<?( x ) ~ <f>(x)<f>h (y,z), ho (y, *) = ^ e 



■••1 



(2.23) 
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Plugging (2.23) into (1.5), multiplying both sides by (j>ho(y,z) an d integrating over 
(y, z) el 2 , we get 



ji 4>(x) = --(f> xx + V]_(x)(f> + Pi\(p\ 2 (p, -oo < x < oo, 



(2.24) 



where 



2 

Vi(x) = ^ fi = pf \<f> ho (y,z)\ 4 dydz = 

2 JR2 



2tt 



Using the results in the Appendix A for ID GPE, again we get approximate ground 
state in three typical regimes: 

Regime I. Weakly interacting regime, i.e. /3° = l3^J^ y ^ z /2Ti = o(l), the ground state 
is approximated by the harmonic oscillator ground state: 

9 (x) « ho (x)0 ho (y, z) = (f) ho (x, y,z), xel 3 , ly > 1 & 7z > 1 & = o(l), (2.25) 



7y + 72 + i + W) = + J + O (/*(7»7.) 1/a ) , 



9 2 2 2 2 

^ - + * + W) = + ^ + O (/?( 7w7 ,) 1/2 ) • 



(2.26) 
(2.27) 



Regime II. Intermediate interacting regime, i.e. /3" 
ground state can be approximated by 

9 (x)^0f(x) :=0j D (z)0 ho (y,z), 



Eg ~ 

where 



= PVW~*/ 27T = 0(1), the 
x G R 3 , (2.28) 

(2.29) 



-4- '"V ~ 



2 ™ 9 7 2 

^ CS == M0i D (*)<My, z)) = ?^p± + ^ 1d(0 ^) := 2V±^i + ^ (2 30) 



ID 



E 



4" 



roo 


1 




J — oo 


2 




roo 


1 


^ D (x) 


I —oo 


2 


dx 



+ V l {x)\ (j) f{x)\* + ^\<pf{x)\* 



+ V 1 (x)\^(x)\ 2 + Pt\C^) 



dx, 
dx. 



Here 0* D , and are the ground state, energy and chemical potential of the 
ID problem (2.24). In this case, one only needs to solve a ID problem numerically 
and thus significantly save computational time, memory and cost. 

To verify (2.13), (2.14) and (2.15) numericallly, Table 4 lists the error ||0 s (x) - 

(//? s (x)|| L 2, and Figure 3 shows the error ^ 9 ^ g - , for different (3 and 7 := -y y = 7z . 
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1/7 

/ / 




1/12.5 


1/25 


1/50 


1/100 


1/200 


3 = 


25 


0.1512 


0.1283 


0.1076 


0895 


0.074 


rate 






0.24 


0.25 


0.27 


0.28 


3 — 


50 


2232 


1914 


0.1162 


1363 


1 1 36 


rate 






0.22 


0.24 


0.25 


0.26 


P = 


100 


0.3150 


0.2742 


0.2357 


0.2006 


0.1692 


rate 






0.20 


0.22 


0.23 


0.25 


P = 


200 


0.4228 


0.3740 


0.3269 


0.2826 


0.2418 


rate 






0.18 


0.19 


0.21 


0.22 


P = 


400 


0.5389 


0.4851 


0.4316 


0.3798 


0.3309 


rate 






0.15 


0.17 


0.18 


0.20 



Table 4: Error analysis of \\<f> g - 
trap. 



cCSl 



L 2 for the ground state in 3D with a cigar-shaped 





-e- 


Y=12.5 




Y=25 




y=50 




y=100 




Y=200 



(b) 



3.5 



4.5 

ln(P) 



5.5 



\E — E cs \ 

Figure 3: Convergence rate of E 9 in 3D with a cigar-shaped trap: (a) With 
respect to 7; (b) with respect to (3. 



From Tab. 4, Fig. 3 and additional numerical results in [18], when 7 = 7^ = 
7 Z ^> 1, ft" — P r y/2'K = 0(1) or 1, and /?7 _1 = o(l), we can draw the following 
conclusion: 
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I csi n (o 1/3, \ n (^!U^l\ 
k " I = O ((3 7 1 ln 7 J , = O y ^ 2/3 j . 

On the contrary, when 7 = 7^ = 7^ >> 1, /3" = j3 r y/2ii = 0(1) or ^> 1, and 
/9 7 _1 = 0(1) or ^> 1, the errors do not decrease when 7 increases. This suggests 
that one cannot reduce the 3D GPE (1.5) to ID GPE (2.24) when f3 > 1 and 7 > 1 
with p-f" 1 = 0(1) or > 1. 

Regime III. Strong interacting regime, /3" = (3^y^%/2-K 3> 1, noticing (A. 17) with 
d = 1, the ground state in 3D is approximated by the multiplication of the TF 
approximation in ^-direction and the harmonic oscillator approximation in yz-plane: 



0,(x) » 0j F2 (x) := (/^(x^y, z): x G K 3 , (2.31) 



where 

tp (x) = { ^ - *72) - 2 < 2/,? D F , (2 32) 

[ otherwise, 

;;TF _ 1 /3^\ 2/3 _ m 2/3 (lyl Z ) 1/3 (? ^ 

^ m ~2\Y) ~ 2(4^ • 

Plugging (2.28), (2.24), (A.7) with d = 1, (2.33), (A.16) with d = 1 and (3 X = (3<t 
into (2.14), we get the approximate energy: 

E g = E{ct )g )=E a { ( t ) f{x)^Uy^)) + 0(h l jHn ly ) 

-H^-H (f )" + Onff + CO + O ln 7 , 

« ^ TF2 + 0(/? 7 ; /3 ln7,), (2.34) 

E ° ~^2~ + 10(4vr)2/3 • (2 ' 35) 
Similarly, we get the approximate chemical potential: 

fi g ^^ F2 + 0(f3^ 3 \n ly ), (2.36) 

where 



where 



TF2 _ Jy-riz o \h> lyjz) (r> o 7 \ 



TF2 _ 7, + 7* , 3 2 / 3 (/3 2 7 ,7 2 ) 1/3 
2 2(4tt) 2 / 

Specifically, if 7^ = 7 2 := 7, then (2.33), (2.34) and (2.36) collapse to 

^< F2 + 0(/37 1/3 ln 7 ), /i^/iJ F2 + 0(/57 1/3 ln7), (2.38) 

^ 1D " 2(4vr)2/3' ^ " 7+ l (4vr)2/3 ' ^ " 7+ 2 (4vr)2/3 " 
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1/7, 




1/12.5 


1/25 


1/50 


1/100 


1/400 


8 = 


25 


2360 


0.1958 


0.1606 


0.1309 


0.1061 


rate 






0.27 


0.29 


0.30 


0.30 


8 — 


50 


3385 


2851 


2373 


1958 


1604 


rate 






0.25 


0.26 


0.28 


0.29 


P = 


100 


0.4691 


0.4022 


0.3407 


0.2857 


0.2375 


rate 






0.22 


0.24 


0.25 


0.27 


P = 


200 


0.6212 


0.5440 


0.4706 


0.4026 


0.3408 


rate 






0.19 


0.21 


0.23 


0.24 


P = 


400 


0.7856 


0.7031 


0.6221 


0.5442 


0.4707 


rate 






0.16 


0.18 


0.19 


0.21 



Table 5: Error analysis of ||(/> 2 
cigar-shaped trap. 



! (x)) 2 || L i for the ground state in 3D with a 



1/7, 


1/12.5 


1/25 


1/50 


1/100 


1/200 


(3 = 25 


0.1226 


0.1447 


0.1711 


0.2025 


0.2390 


P = 50 


0.1633 


0.1985 


0.2402 


0.2880 


0.3426 


p = 100 


0.2110 


0.2597 


0.3158 


0.3794 


0.4516 


P = 200 


0.2517 


0.3100 


0.3769 


0.4531 


0.5393 


P = 400 


0.2772 


0.3410 


0.4146 


0.5001 


0.5975 



Table 6: Error analysis of \\(f> g 
shaped trap. 



lTF2\ 



L oo for the ground state in 3D with a disk- 



To verify (2.31), (2.38) and (2.39) numericallly, Tables 5 and 6 list the errors 
||0 2 (x) — (0j F2 (x)) 2 || L i and ||0 2 (x) — (0j F2 (x)) 2 ||x,oo, respectively, and Figure 4 shows 

\E — _B TF2 | 

the error „ 3 — -, for different P and 7. 



E a 



From Tabs. 5-6, Fig. 4 and additional numerical results in [18], when 7 := 7 y - 
7 2 ^> 1, Pi — P r y/2n ^> 1 and fij^ 1 = o(l), we can draw the following conclusion: 



II*,- 



L 2 



o 



'p 1 / 3 In 



7 



7 



1/3 



/TF2\2|| 



o 



'p 1 ! 3 In 



7 



7 



1/3 



\E n - E. 



TF2i 



\Hg~H 



TF2i 

9 



~-0(p 7 1/3 In 7 
0(p 1 1 ' 3 \n 1 ), 



\E n 



E TF2 \ 



E 9 

\»9 -/^J F2 1 



o 



'P 1 ' 3 In 



7 



7 



2/3 



= O 



'P 1 / 3 ln 7 N 



7 



2/3 



These results confirm the asymptotic results (2.38), (2.39), (2.34) and (2.36). Fur- 
thermore, the additional numerical results in [7, 18] indicate that 0^ F2 (x) does not 
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converge pointwise to the ground state ff (x) when j z — > oo and (5 > (cf. Tab. 6). 

2.4 Choice of initial data for computing ground states of 3D 
GPE 

In order to numerically compute the ground state of the 3D GPE (1.5) efficiently, 
for a given numerical method, an appropriate choice of the initial data is also very 
important such that the iterative number is highly reduced and thus the total com- 
putational cost is greatly saved [7, 4, 13]. Based on the discussion in the previous 
subsections, here we provide an appropriate choice for the initial data for different 
parameter regimes f3, 7 y and ^ z in (1.5): 

• I. When i) i y = 0{ l), 7, = 0(1) and (3 = o(l) or 0(1); ii) ly = 0(1), lz > 1 
and ^ = (3^ 1z /2tt = o(l); or iii) ly > 1, lz > 1 and /3« = f3^r z /2n = o(l); 
an appropriate choice of the initial data is 

0o(x) = Mx,y,z) = U y 3/ 4 e 5 > x G R ■ 

• II. When ly = 0(1), lz > 1, /3 2 a = ^ lz /2n = 0(1) and /? 7 ," 3/2 = o(l), an 
appropriate choice of the initial data is 

0o(x) = 0f ho (^) = ^ e- 2 ^ 0f(x,y), x G K 3 ; 

where (j)g D (x, y) is the ground state of the 2D GPE (2.9) and can be computed 
numerically [7, 4, 13]. 
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III. When ly = 0(1), 7z > 1, ^ = (3^ z /2n > 1 and (3% z l 2 = o(l), an 
appropriate choice of the initial data is 

o (x) = (fiE(x, y) <j) ho (z) = -f^ e" V- y), x G K 3 ; 

where y) is given in (2.17). 

IV. When 7 , » 1, 7z » 1, /3? = P^%T z /2n = 0(1) and /5( 7y 7,)- 1/2 = o(l), 
an appropriate choice of the initial data is 



o (x) = 4>f{x) ho (y, z) = { ^^- e" 2 ^ 2 ^ <C(x), x G K 3 ; 



where <p l ^(x) is the ground state of the ID GPE (2.24) and can be computed 
numerically [7, 4, 13]. 

V. When ly > 1, lz > 1, = f3^yul^ > 1 and /3(7 y 7z)~ 1/2 = o(l), an 
appropriate choice of the initial data is 



O (X) = <f>™(z) <f> ho (y, Z) = e -Z^±Z-! TF (x)) x G R3; 



where (ff^(x) is given in (2.32). 

• VI. For all other good choice of the initial data is 

0o(x)=0j F (x), xGK 3 ; 

where 0j F (x) is given in (2.4). In fact, in these cases, the 3D GPE (1.5) cannot 
be reduced to lower dimensions! 

Our numerical experiments showed that the initial data in cases II- V are much 
better than initial data often used in the current literatures. 



3 Dimension reduction for time-dependent GPE 

In this section, we will discuss dimension reduction of the 3D time-dependent GPE 
(1.3) in certain limiting frequency regimes and provide convergence rate of the re- 
duction based on our numerical results. 

3.1 In disk-shaped condensates 

In the disk-shaped condensates, i.e. for uj x ~ u> y and u z ^> uj x (■<=>• 7 y ~ 1 and 
7 2 3> 1), the 3D GPE (1.3) can be reduced to 2D GPE by assuming that the time 
evolution does not cause excitations along the z-axis since they have a larger energy 
of at least approximately hu z compared to excitations along the x- and y-axis with 
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1/7* 




1/25 


1/100 


1/400 


1/1600 


3 = 


1 


1 954E-3 


8 619E-4 


3.547E-4 


1.346E-4 


rate 






0.59 


0.64 


0.70 


8 — 


If) 


1 056E-2 


3 968E-3 


1.446E-3 


5.149E-4 


rate 






0.71 


0.73 


0.74 


P = 


100 


3.709E-2 


1.339E-2 


4.769E-3 


1.689E-3 


rate 






0.74 


0.75 


0.75 


P = 


1000 


1.132E-1 


4.216E-2 


1.511E-2 


5.364E-3 


rate 






0.71 


0.74 


0.75 


P = 


10000 


2.902E-1 


1.256E-1 


4.707E-2 


1.687E-2 


rate 






0.60 


0.71 


0.74 



Table 7: Error analysis of ||03 — 0ho||i 2 for dimension reduction from 3D to 2D. 

energies of about hu x . Following the physics literature [23, 17, 7, 5, 21], for any 
fixed (3 > and when 7 Z — > oo, we assume that the condensation wave function 
along the z-axis is always well described by the ground state wave function which is 
well approximated by the harmonic oscillator in z-direction and set [21, 17, 7, 5] 



if> = ip 2 (x,y,t)(j) 3 (z), (j) 3 (z) = |0 9 (x)| 2 dxdy^j m (j) ho (z) 

(3.1) 

where <j) g (x,y,z) defined in (1.9) is the ground state solution of the 3D GPE (1.3). 
Following the procedure used in [21, 17], the 3D GPE (1.3) can be reduced to a 2D 
GPE with x = (x,y): 

1 = -\ V2 ^^ + ^(*,y)iKx,t) + ^|^(x,t)|V(x,t), (3.2) 

where 

fo = P J™jt(z)dz^3 a 2 -.= B^ v 2 (x,y) = ^(x 2 + 1 y). (3.3) 

To verify (3.1) and (3.3) numerically, Table 7 lists the error ||0 3 (^) — (f)ho(z)\\ L 2, 
and Figure 5 shows the error ^ 2 ^ 2 vs. ^ z and (3, for different (3 and 7 Z . 

From Tab. 7, Fig. 5 and additional numerical results in [18], for fixed 3 > and 
when 7 Z — > oo, we can draw the following conclusion: 



\\H*)-<h»(*)\\* = o(^' 1 ' 



z 



3/4 

7* 



Ik- = O f ^M" 7 ' ) > 11^) -^oWIUi = O 



7* / \ lz 
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Figure 5: Convergence rate of ^ 2 for dimension reduction from 3D to 2D: (a) 
With respect to r y z ; (b) with respect to (5. 



In addition, we compare the errors between the solutions of the 3D GPE (1.3) 
and its 2D reduction (3.2) for different 7 2 . In order to do so, for any given j z , let 
vjj 3D (x,y, z,t) be the numerical solution of the 3D GPE (1.3) with ^ x = % = 2, 
P = 1 and the initial data i/)(x, y, z, 0) = ipo(x,y,z) chosen as the ground state of 
(1.3) with 7 X = 7 y = 1, /3 — 1 [4, 7]. This 3D dynamics corresponds to a BEC 
which is in its ground state when at t — the trap frequencies in x and y direction 
are suddently doubled, i.e. setting = l y = 2. Similarly, let ip 2B (x,y,t) be the 
numerical solution of the 2D GPE (3.2) with ^ x = 2, ^ y = 2, (3 2 = (3% := and 
initial data ip(x, y, 0) = ipo(x, y) chosen as the ground state of (3.2) with j x = j y = 1, 
@2 = (3% = Pyf%H [4, 7]. Again, this 2D dynamics corresponds to a BEC which is 
in its ground state when at t — the trap frequencies in x and y direction are 
suddently doubled, i.e. setting 'jx = ly = 2. In fact, ip 2D is the solution of the 2D 
reduction problem. In order to do the comparison, we introduce 



03(M) = 

3D (x,t) 



vjj 3B (x,y, z,t) dxdy 



1/2 



0ho(^) 



7. 



1/4 



x,£):=-0 (x,y,t)(f) ho (z), 



x G 



71-1/4 
p3 



7 z z 
2 



(3.4) 
(3.5) 



and the condensate widths 



crJt) = f a 2 ^ 3D (x,t) 2 <ix, cr"(i) = / a 2 ^ DS (x, t) 2 dx, a = x,y,z 

JM 3 JM 3 



(3.6) 

Figure 6 shows the errors HV^z, — <l>h.o(z)\\L°°, — °""| = Wy — <r%\, &z — c z = 
|^ 3D (0, t)| 2 — |-?/> DS (0, t)\ 2 for different 7 Z . Here the numerical solution 
ifj 3B and ip 2D are obtained by the time splitting spectral method [5, 3]. 



a z — j and 
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From Fig. 6, the dimension reduction from 3D time- dependent GPE (1.3) to 2D 
GPE (3.2) when f3 = 0(1) and 7 Z 3> 1 is verified numerically. Furthermore, we have 
the following convergence rate: 

\\Mz,t) ~ = O a x (t) = a a x (t) + (-^J , lz > 1, 

= i + °(^74) > I^CM)! 2 = l^ DS (x,t)| 2 + ^ . 
3.2 In cigar-shaped condensates 

Similarly, in the cigar-shaped condensates, i.e. uj v ^> uj x and oj z ^> u x ( < ^=>- 7 y ^> 1 
and 7 Z ^> 1), for any fixed f3 > and when 7 y — > oo and 7^ — > oo, the 3D GPE 
(1.3) can be reduced to ID GPE with x = x [21, 23, 5, 7, 17]: 

. ^t)_ = _l dxx ^ Xit) + v^xMxrf+faM^tWfat), (3.7) 

where 



& = |^> 23 (2/, ^)| 4 ^2/^ ^ := V 1 (x) = Z-, (3.8) 
23 (y,z) = (j^ |0,(x,y,z)| 2 cfej » ho (^) := Uy ^ /2 e 5—. (3.9) 



Similarly, to verify (3.8) and (3.9) numerically with 7 := 7 y = 7z, Table 8 lists 
the error lMg^Mlz^h°fe)lk2S an( j Figure 7 shows the errors ^ ^ vs. 7 and /?, for 
different 7 and (3. 

From Tab. 8, Fig. 7 and additional numerical results in [18], for any fixed 
(3 = 0(1) and when 7 := 7^ = j z — > 00, we can draw the following conclusion: 

||/2 / \ ,2 M || n^V3 n/ l/3 lTln A H023(-)-^ho(-)lUi n ^ 1/3 ln 7 ^ 

ll^s(0 - ^ho(OIUi = O [p 7 ln 7 J , 1|^3(.)|| £1 = \—fJ^ ) ■ 

3.3 For vortex interaction in 3D 

In this subsection, we numerically study dimension reduction for vortex interaction 
in 3D disk-shaped condensates. In this case, the assumption (ps(z) ss <f>h (z) in (3.1) 
is no longer valid. In order to do so, for any given 7^, let 0f (x, y, z) be the central 
vortex line state with winding number m — 1 of the 3D GPE (1.3) with 7^ = 7 y = 1 
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1/7 

/ / 




1/12.5 


1/25 


1/50 


1/100 


1/400 


3 = 


25 


0.1727 


0.1412 


0.1145 


0.0924 


0.0743 


rate 






0.29 


0.30 


0.31 


0.32 


3 — 


50 


2591 


2135 


1 746 


1419 


1 148 

W . J. 1t:U 


rate 






0.28 


0.29 


0.30 


0.31 


P = 


100 


0.3378 


0.3151 


0.2606 


0.2141 


0.1748 


rate 






0.26 


0.28 


0.29 


0.29 


P = 


200 


0.5334 


0.4517 


0.3791 


0.3156 


0.2608 


rate 






0.24 


0.25 


0.26 


0.28 


P = 


400 


0.7285 


0.6266 


0.5345 


0.4521 


0.3792 


rate 






0.22 


0.23 


0.24 


0.25 




(a) - |n <^> (b) ln <e> 

Figure 7: Convergence rate of for dimension reduction from 3D to ID: (a) 

With respect to 7 = 7^ = 7 Z ; (b) with respect to j3. 



and (5 — 10 [4, 8]; and ip 3D (x,y, z,t) be the numerical solution of the 3D GPE (1.3) 
with 7 X = 7^ = 1 and (5 — 10 and the initial data if)(x, y, z, 0) = ipQ D (x, y, z) chosen 
as 

,3D/ ^_ <j>f {x-x ,y,z)<j)l D {x + x ,y,z) 3 

ro \ x iVi z ) — 11 / 3D / \ J.3D/ 1 Ml ' \ x iVi z ) fc M ) 

W\\ x - x , y, z)(j)'l u {x + x , y, z) \ L i 

where x > is a constant. This 3D dynamics of BEC corresponds to the interaction 
of two vortex lines in 3D. Formally, when 7 2 3> 1 [4, 8, 5], 

ft \ 1 / 2 (2>y W 4 

M*) = [J R2 \^l D (x,y,z)\ 2 dxdyj » e~^ z ± ho (z). 
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Similarly, let 4>1 D (x, y, z) be the central vortex state with winding number m — 1 of 
the 2D GPE (3.2) with lx = ly = I and (5 2 = [4, 8]; and ip 2D {x,y,t) be the 

numerical solution of the 2D GPE (3.2) with 7x = 7y = 1 and /3 2 — P\f^ an d the 
initial data ip(x,y,0) = ^^{x^y) chosen as 



i2D 



x 



- x ,y)(f) 2D (x + x 0: y) 



Again, this 2D dynamics of BEC corresponds to the interaction of two vortices in 2D. 
Let x 2 d(£) = (x(t),y(t)) T be the location at time t of center of the vortex initially 
located at (xq,0) in the 2D dynamics, and x 3 d(£) = (x(t),y(t)) T be the location in 
the xy-pl&ne at time t of center of the vortex initially located at (xq, 0, 0) in the 3D 
dynamics, i.e. center of the vortex from the wave function tp 3B (x,y,0,t). Figure 8 
shows the trajectory of x 2D (t) and x 3D (i) with x — 1 for different 7 2 . 

From Fig. 8, for fixed 7 2 , the trajectories of the vortex centers in 2D and 3D 
dynamics agree qualitatively (cf. Fig. 8). But for fixed j3 = 0(1), when 7 2 — > oo, 
the larger is 7 Z , the larger is the error. This is because for larger 7 Z the condensate 
becomes flat and thus the vortex line dynamics in 3D bent more frequently which 
causes the oscillatory nature in the trajectory. 

For dimension reduction of the GPE with general initial data, we refer to [6, 10, 

30]. 



4 Conclusion 

Dimension reduction of the three-dimensional (3D) Gross-Pitaevskii equation (GPE) 
for Bose-Einstein condensation under different limiting parameter regimes was stud- 
ied asymptotically and numerically. For ground state of the 3D GPE, we found the 
convergence rate of the reduction and provided the approximate energy and chemi- 
cal potential in both disk-shaped condensation and cigar-shaped condensation. Our 
extensive numerical results confirmed the reduction and the convergence rate. In 
addition, we identified the parameter regimes in which the reduction is invalid. For 
dynamics of the 3D GPE, our numerical results confirmed the reduction and pro- 
vided convergence rates in certain limiting parameter regime. 



Appendix: Energy and chemical potential approximations for the 

ground states 

The 3D time-independent GPE (1.5), 2D GPE (2.9) and ID GPE (2.24) can be 
written in a unified way [17, 5, 7] 

H 0(x) = -^V 2 0(x) + V^(x)0(x) + /3,|0(x)| 2 0(x), x G R d , (A.l) 
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>> 




(b) 










(c) x ( d ) 

Figure 8: Trajectory of the vortex centers in 2D and 3D dynamics. Blue line (or 
flat curve): X2d(£); and red line (or oscillatory curve): X3d(£). 



under the normalization condition 

|2 . 



fx)| 2 dx = 1; 



(A.2) 



where (3% — (3, fa = /3f , Pi = and V^x) = V^(x). The energy functional is defined 

as 



E(cf>) 



i|V0(x)| 2 + K(x)|0(x)| 2 + ^| 



dx. 



(A.3) 



Al. Thomas-Fermi (TF) approximate energy and chemical potential. 

When (3 d > 1, and -y y = 0(1) with d = 2 or j y = 0(1) and -f z = 0(1) with d = 3 in 
(A.l), we can ignore the kinetic term and derive the TF approximation: 



m7 07 (x) 



^(x)0f(x)+/3 d |0f(x)| 2 0f(x), 



x G 



(A.4) 
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Solving (A. 4), we obtain the TF approximation for the ground state: 



C(x) = ^ - v«*)) /A, *«*) < „J F , (A 5) 

[ 0, otherwise. 
Plugging (A. 5) into (A. 2) with = FF , we obtain [5, 7] 



TF 1 



^) 2/3 , 1=1, 

1/2 



2 



f^T , d = 2, (A.6) 



/ 15/3 7m 7z \ 2/5 



d = 3. 



Since 0j F (x) is not differentiable at Vd(x) = yuJ F , as observed in [7, 5, 8], _E(0j F ) = 
oo, thus one cannot use the definition (A. 3) to define the energy of the TF approx- 
imation (A. 5). Therefore, noticing (1.6) and (1.8), as observed in [8], here we used 
the way to calculate it: 

ej f * Eg = E {<t> g ) = M0 9 ) - £in t (<y - »7 - Euf 9 F ) = j^7- ( A - 7 ) 

A2. First-order approximate energy and chemical potential. 

When 7 y = 1 with d = 2 or ^ y = 7 2 = 1 with d = 3 in (1.5), the ground state of 
the nonlinear eigenvalue problem (A.l) is symmetric, i.e. 9 (x) = 0(r) with r = |x|, 
and satisfies: 

1 t\ f rd ^d^r)\ + {Vd{r) _ ^ 0(r) + ^ 03(r) = Q) o < r < oo, (A.8) 



2r d 1 dr \ dr 

where Vd(r) = r 2 /2. Following the method used in [29] for d — 3, we choose 
R = y / 2 / uJ F such that V d {R) = f = /iJ F When |r - R\ < i?, we have 

W - n « W - „J F = £ - ^ = (r - i?)^ « (r - R)R. (A.9) 

Noticing (A.9) and dropping the first order term in (A.8), we obtain 

1 d 2 <j)(r) 



2 dr 2 



(r-R)R(j)(r)+(3 d (f) 3 (r)=0, < r < oo. (A.10) 



Introducing a change of variables, s = (r — R)/l and 0(r) = a<f>(s) with 2i?/ 3 = 1 
and 2(3 d a 2 l 2 = 1, we can reduce (A. 10) to 

4>"(s) - (s + 2 (s))0(s) = 0, -oo<s<oo. (A.ll) 

For the solution of (A.ll), as s — > +oo, — > 0, dropping 3 term in (A.ll), we have 

0"(s) - s4>(s) = 0. (A.12) 
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Thus we have the asymptotics for the solution 



4>(s — > +00) 



A 



-S.s 2 /3 



2sV4 



e 3- 



0.794. 



On the other hand, as s — > —00, dropping 0"(s) in (A. 11), we have 

s + 2 (s) = 0. 

Thus we get 



(A.13) 
(A.14) 

4>(s -> -00) ~ v^s- (A.15) 
Choosing e such that I <^ e <^ R, using 9 ps 0j F for r e [0, it! — e] and 
9 ~ a0((r — R)/l) for r e [R — e,oo), plugging (A. 5), (A.15) and (A. 6) with 
1y = lz = 1 into (1.8), we get the approximate kinetic energy of the ground state 

pR—e roo 

/ (0>)) 2 r d_1 dr + / (<p'Jr)) 2 r^ 1 dr 

JO JR—e 

'l 
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Cd 
2 



W\ 2 



dr 



/c 



2 



1 r R ~ e 



„d+l 



2/3, Vo 2/iJF - r 2 



dr + 



J or 



2 

0* 
2 

Cd 
2 

Cd 
2 

cy^ 



1 

I 

R d 



R-e 



2(3 d Jo 2/iTF _ r 2 



ln^- + D d ) + 



dr + 



/oo 

/oo 



<f>'(s) (ls + R) d 'Hds 

d-1 



I J-e/l 
e/l 

a 2 R d ' 1 roo 



ds 



ds 



R d f, R - , 

R d 



I J-e/l 

a 2 R d - 1 



/oo _ 2 / \ 

{ (f)'(s) v / TT^ 2 d In (s + vTT^ 2 ) 



4/ 



lnf + C 



In (2 1 / 3 i? 4 / 3 ) + D d + C] 



,2/(d+2) 
d 



(ln^ + G,), 



(A.16) 



where 



C = -4 f° ln(v / TT^ 2 + s)^-fv / rT7 2 "(0 / ( s )) 

J-e/i v 7 as L 



In (vT 



s 2 + s 



(is 



(is 



VI + s 2 (0'(s)) 2 ds^ 0.706, 



i? 



2/iJF 



((d+l) 2 -l )(3 d 

a. 



1 V(d+2) 
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C^\ {d+ If _1)W2) 

Q " 6(dT2) ' rf " 1 ' 2 ' 3 ' 

Gd = ln^±itl + ^(ln2 + 3^ + 3C), 
W 4 



f 2, 

2tt, 

k 47T, 



' ln4-2, 

-1, 
k In 4 -8/3, 



d = l 
d = 2 
d = 3. 



From (A. 3), (A. 7), (1.6), (1.8), (A. 4) and (A. 16), we can get the first order approx- 
imation for Eg and ji g when (3d >> 1: 



-ttf ■ - - d -\- 2, 



2(d + A) 



{{d + lf-l)(3 d 



C d 



2/(d+2) 



+ 



^2/(d+2) 



rf + 2 

2(rf + 4) 



((rf+l) 2 -l)& 
C r , 



2/(d+2) 



+ 



In 



2/(d+2) I 



(ln/3 d + G d ) 



(A.17) 



/iJ F + E kin (0 9 



l) 2 - 



1 2/(d+2) 



((d+l) 2 -l)& 



2/(d+2) 



+ 



ln& 



4 



2/(d+2) 
d 



,2/"(d+2) (^^d + Cd) 

(A.18) 



These asymptotic results were confirmed by the numerical results in [7]. 
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